function coherency_phase_spa_13(wave_phase,wave_coher,period,time_out,coi_out,lat,time,ilat,itps,anom,anomOC)

close all
figure;
%set(gcf,'Visible','off')          

Yticks = 2.^(0:fix(log(max(period))/log(2)));
load map_phase
colormap(map_phase)



% --------------------------------
subplot(2,2,1:2)

[AX,H1,H2] = plotyy(time_out,anom,time_out,anomOC,'plot');

%You can use the handles returned by plotyy to label the axes and set the line styles used for plotting. With the axes handles you can specify the YLabel properties of the left- and right-side y-axis:
set(get(AX(1),'Ylabel'),'String','SLA Anomaly (cm)') 
set(get(AX(2),'Ylabel'),'String','Chl Anomaly (mgChl.m^{-3})') 
set(AX(1),'YColor','r')
set(AX(2),'YColor','b')
set(AX(1),'Xlim',([time_out(1) max(time_out)]))
set(AX(2),'Xlim',([time_out(1) max(time_out)]))
set(H1,'Color','r')
set(H2,'Color','b')
set(AX,'XColor','k')
xlabel('Distance (km)')
title('a) SLA et Chlorophylle - 30 sud - Septembre 1999'); 
hold off


subplot(2,2,3)
pcolor(time_out,period,wave_phase')
caxis([-pi pi])
shading flat
set(gca,'YScale','log',...
    'YDir','reverse',...
    'YTick',Yticks,...
    'Yticklabel',Yticks,...
    'Tickdir','out');
fig=colorbar('vert');
ylabel('Wavelength (km)')
title('b) Phases');
hold on
plot(time_out,coi_out,'k')   % cone-of-influence, anything "below" is dubious
hold off            
xlabel('Longitude');

subplot(2,2,4)
pcolor(time_out,period,wave_coher')
shading flat
set(gca,'YScale','log',...
    'YDir','reverse',...
    'YTick',Yticks,...
    'Yticklabel',Yticks,...
    'Tickdir','out');
fig=colorbar('vert');
%ylabel('Longueurs d''onde (km)')
title('c) Coherency');
hold on
plot(time_out,coi_out,'k')   % cone-of-influence, anything "below" is dubious
hold off
xlabel('Distance');



%cd /data/couleur/travail_en_cours/ELODIE/RW_ANALYSE/FIG_CROSS/
%print('-f1','-depsc2',['./COHER_spa_SLA_OC_crlog_wwfilt_',num2str(lat(ilat)),'_',num2str(time(itps)),'.eps'])
%close all
%cd /home/couleur/gutknech/ANALYSE_RW/METRICS_RW/FIG_CROSS/      
